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ABSTRACT 


This  thesis  presents  a  novel  antenna  structure  that  satisfies  the  challenging 
requirements  of  an  air  coupled  high  speed  ground  penetrating  radar  (GPR).  The 
desired  GPR  system  is  to  achieve  high  spatial  resolution  and  accurate  inspection 
results  while  scanning  at  relatively  high  speed  for  highway  pavement  and  bridge 
deck  inspection.  This  work  utilizes  the  Ultra  Wide  Band  (UWB)  antenna  design  to 
achieve  both  physical  and  electrical  requirements  imposed. 

The  design  procedure  starts  with  a  short  survey  to  discuss  typical  UWB 
antennas  used  for  GPR  applications,  and  various  tradeoffs  of  each  type  specifically 
when  used  for  Air  Coupled  GPR  applications.  Our  structure  anatomy  is  presented, 
followed  by  a  theory  introduction  that  mainly  focuses  on  achieving  good  impedance 
matching  throughout  the  proposed  antenna  structure.  A  proof-of-concept  MATLAB 
model  is  created  to  evaluate  the  preliminary  physical  dimensions  that  can  achieve 
minimum  reflections  at  antenna’s  feed  point.  These  dimensions  are  then  used  in 
SolidWorks  to  create  a  3D  model  that  is  imported  later  in  HFSS  to  obtain  accurate 
electromagnetic  characteristics.  Furthermore,  fine  tunings  are  performed  to  the 
antenna  structure  to  optimize  both  gain  and  impedance  matching.  The  SolidWorks 
3-D  structural  model  is  finally  used  for  antenna  fabrication.  The  measurements 
recorded  from  the  field  experiments  using  the  prototypes  manufactured  are 
compared  to  the  simulation  results  confirming  our  initial  findings.  Both 
measurements  and  simulation  results  demonstrated  very  small  reflection  loss  across 
the  700  MHz  ~  6  GHz  frequency  band  with  a  very  high  directed  gain  and  radiation 
efficiency. 
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CHAPTER  1:  INTRODUCTION 


The  primary  objective  of  this  thesis  is  to  focus  on  one  of  the  essential  blocks  in 
any  Ground  Penetrating  Radar  (GPR)  systems.  We  set  our  target  on  coming  up  with  a 
systematic  approach  for  antenna  design,  modeling,  manufacturing,  and  field  testing.  A 
good  base  was  already  established  earlier  by  my  colleagues  at  the  University  of 
Vermont  to  investigate  the  feasibility  of  an  Air-Coupled  GPR  system  that  can  be  used 
in  bridge  deck  or  highway  subsurface  inspections  without  disrupting  traffic. 

In  the  year  2000,  Jing  Qiong  Hu  led  the  way  by  first  introducing  the  Good 
Impedance  Matching  (GIMA)  antenna,  which  demonstrated  low  reflections  at 
aperture  and  proved  the  theoretical  and  practical  feasibility  of  using  GIMA  antennas 
in  a  GPR  module.  Later  on,  in  2012,  Xianlei  Xu  et  al.  introduced  the  dual-band  dual¬ 
channel  air-launched  UWB  GPR  to  facilitate  highway  and  bridge  inspection.  The  new 
novel  system  achieves  high  pulse  amplitude  and  low-level  ringing  that  also  improves 
the  survey  speed  by  using  the  8-GSPS  real-time  data  acquisition  unit  to  accelerate 
reflection  sampling.  The  promising  results  assure  that  this  system  can  achieve  good 
resolution  while  scanning  on  highway  driving  speeds.  However,  the  new  system 
revealed  the  need  for  a  new  antenna  design  that  can  achieve  higher  gain  at  a 
frequency  band  500MHz  to  6GHz  with  reflections  at  feed  point  Sll  less  than  -lOdB. 
And  given  that  the  system  will  ultimately  be  installed  under  a  truck,  we  have  to 
understand  the  effect  of  the  vertical  length  of  the  antenna  on  the  overall  performance. 
These  factors,  in  addition  to  design  mechanical  reliability  for  long  term  operations, 
exposed  the  need  to  develop  a  systematic  approach  to  determine  the  critical 
dimensions  that  contributes  to  the  antenna  design  and  model  the  effect  of  each 
dimensional  variation  on  the  targeted  performance. 
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In  this  chapter  we  will  cover  the  required  literature  survey  that  will  first 
introduce  the  GPR  systems:  their  typical  applications,  what  makes  up  a  good  GPR 
system,  and  the  different  types  and  trade-offs  in  the  system  design.  Afterwards,  we 
move  to  the  antenna  introduction  including:  the  ultra-wide  band  antenna,  why  it  is 
suitable  for  ground  penetrating  radars,  and  what  types  of  UWB  antennas  are  suitable 
for  an  Air-Coupled  GPR.  And  finally  we  review  the  system  our  antenna  will  take  part 
in  and  explore  the  system  requirements  imposed  on  the  antenna  structure. 

1.1.  GPR  literature  survey 

The  Ground  Penetrating  Radar  (GPR)  technology  is  increasingly  being  used  in 
non-destructive  testing  and  through-wall  imaging.  Various  commercial,  governmental 
and  military  applications  often  require  the  need  to  detect  buried  objects  or  reveal  their 
nature.  Sometimes  it  is  also  required  to  determine  the  internal  structure  of  materials, 
similar  to  road  layer  thickness,  bridge  construction  and  tunnel  linings.  GPR’s  are  very 
useful  for  investigating  layers’  boundaries  and  can  easily  deliver  onsite  results. 

1.1.1.  Theory  of  operation 

A  typical  surface  image  is  captured  through  transmitting  a  short  impulse,  less 
than  InSec  in  width,  of  electromagnetic  signal  into  the  structure  under  test.  The 
energy  reflected  from  discontinuities  in  material,  which  has  different  electrical 
impedances,  is  collected  by  a  receive  antenna,  processed  by  means  of  signal 
processing,  and  finally  presented  on  a  display  unit.  Figure  1  shows  a  block  diagram  of 
the  system. 
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Figure  1:  Block  diagram  of  a  typical  GPR  system 


In  Figure  2,  the  vertical  depth  of  a  single  point  Pi  located  at  the  boundary  of  two 
materials  having  different  dielectric  constants  can  be  calculated  with  the  aid  of  the 
electromagnetic  wave  propagation  velocity  equation.  Equation  (1)  introduces  the 
propagation  velocity  formula  assuming  a  homogeneous,  isotropic  material,  where  c  is 
the  speed  of  light  or  the  propagation  velocity  in  free  space  and  sri  is  relative 
permittivity  of  the  first  material. 


c 


(1) 


V^rl 


The  depth  to  the  target  point  can  then  be  derived  from  equation  (2),  where  d  is 
the  vertical  depth  of  the  point  Pi,  v  is  the  propagation  velocity  inside  the  first 
dielectric  material,  62  is  the  refraction  angle  that  can  be  easily  calculated  using  Snell’s 
law  given  the  incident  wave  angle  9 /,  and  Vf  is  the  time  difference  between  receiving 
the  first  reflected  pulse  Rxi  and  the  second  reflected  pulse  RX2.  Figure  2  explains  the 
typical  behavior  of  a  transmitted  impulse  Tx  given  that  so<  sri<  sr2. 


2  2  V^7 


(2) 
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Figure  2:  Incident,  refracted,  and  reflected  waves 

The  thickness  of  the  second  layer  buried  within  the  first  material  can  be  easily 
calculated  if  the  depth  of  point  P2  is  found.  A  cross-sectional  image  of  the  material 
can  be  constructed,  point  by  point,  given  that  the  approximate  dielectric  constants  of 
the  different  layers  and  their  refractive  indices  are  known. 

Even  if  we  don’t  know  the  nature  of  the  buried  layer,  we  can  still  acquire  some 
valuable  information  from  its  reflected  signal.  The  reflection  coefficient  T12  is  the 
measure  of  how  much  energy  is  reflected  from  the  second  material  back  to  the  first 
material.  It  can  be  calculated  using  the  characteristic  impedances  ZQi  and  Z02. 


r12 


Z01  zo2 

Zol  +  Z02 


(3) 
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The  characteristic  impedance  ZG  is  defined  as  the  electric  field  over  the  magnetic 


field  inside  the  material,  which  can  be  calculated  using  equation  (4): 


O  —  J(jl)£ 


(4) 


Where: 

G)  is  the  angular  frequency 

p  is  the  material  magnetic  permeability 

s  is  the  material  electrical  permittivity 

a  is  the  material  electrical  conductivity 


If  the  material  is  dielectric,  a  tends  to  be  very  small  in  value  and  the 
characteristic  impedance  can  then  be  approximated  to  Z0  =  yj [i/s 

In  conclusion,  firing  a  short  impulse  of  electromagnetic  wave  towards  a  surface 
not  only  can  detect  the  depths  and  thicknesses  of  its  buried  layers,  but  also  can 
partially  reveal  their  material  characteristics. 


1.1.2.  GPR  types  and  applications 

Scanning  the  surface  of  the  structure  under  test  with  constant  velocity  in  linear 
path  would  develop  a  cross  section  image  along  the  path  scanned.  Repeating  the 
previous  steps  along  parallel  paths  we  can  theoretically  develop  a  3-D  image  of  the 
structure.  One  of  the  early  successful  attempts  to  apply  this  technique  was  developed 
in  the  UK,  by  EVANS  (1963),  to  measure  the  polar  ice  thickness;  probing  depths  of 
several  kilometers  in  the  Antarctic.  A  more  recent  and  relative  to  our  objective  is  the 
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system  developed  by  Warhus  (1993).  He  developed  a  GPR  imaging  system,  using  a 


slot  line  bowtie  antenna,  to  detect  the  locations  for  a  complex  structure  of 
reinforcement  steel  bars  (REBARS)  in  a  grid  form  buried  in  a  concrete  slab. 

According  to  the  diagram  presented  in  Figure  1,  and  based  on  the  antenna 
proximity  to  the  surface  under  test,  ground  penetrating  radars  can  be  one  of  two 
categories;  Ground-Coupled  or  Air-Coupled  GPR.  The  antennas  of  the  Ground- 
Coupled  GPR  are  in  direct  contact,  or  within  very  close  proximity,  to  the  detection 
surface,  which  results  in  smaller  energy  losses  and  higher  detection  sensitivity. 
However,  the  scanning  speed  has  to  be  low  in  order  not  to  expose  the  antennas  to 
possible  damages  if  the  surface  in  contact  has  a  rough  profile.  On  the  other  hand,  Air- 
Coupled  GPR’s  do  not  need  to  be  in  contact  with  the  surface  under  test.  Therefore, 
they  are  useful  in  rapidly  surveying  long  distances  in  relatively  short  time.  For 
example,  GPR’s  have  been  flown  at  heights  of  400m  in  Synthetic  Aperture  Radar 
(SAR)  mode  and  have  discovered  several  buried  metallic  mines.  They  have  also  been 
used  from  earth- orbiting  satellites  to  image  geological  features  buried  beneath  the 
Sahara  Desert. 

In  general,  the  UWB  GPR  is  an  impulse  radar  system  in  which  the  receiver 
collects  the  reflected  electromagnetic  impulses  for  subsurface  feature  identification. 
According  to  D.J.  Daniels  (1996),  the  essential  requirements  that  Ground  Penetrating 
Radars  should  meet  are: 

-  “Efficient  coupling  of  the  transmitted  electromagnetic  energy  into  the  ground”. 

-  “Adequate  penetration  of  the  energy  through  the  ground,  having  regard  to 
target  depth.” 
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-  “Reception  from  buried  objects  or  other  dielectric  discontinuities  of  a 
sufficiently  large  reflected  signal  for  detection  at  or  above  the  ground  surface.” 

-  “An  adequate  bandwidth  in  the  detected  signal,  having  regard  to  the  desired 
resolution  and  noise  levels.” 

-  “An  adequate  signal-to-clutter  ratio.” 

These  factors  have  an  impact  on  the  energy  transmitted  and  the  attenuation  the 
signal  is  suffering  from  while  traveling  in  in  the  media  under  test  and  consequently 
the  received  signal.  The  quality  and  quantity  of  the  received  signal  is  actually  the  key 
indication  of  a  good  GPR  system. 

1.2.  Introduction  to  UWB 

An  antenna  in  general  can  be  described  as  a  transducer  that  guides 
electromagnetic  signal  trapped  in  a  transmission  line  to  radiate  its  electromagnetic 
energy  in  air  or  free  space.  Antennas  may  also  be  considered  as  impedance 
transformers,  that  provides  gradual  coupling  between  input  or  feed  line  impedance  to 
aperture  or  free  space  impedance. 

1.2.1.  Brief  history 

The  fact  is  there  is  no  difference  between  narrow  band  and  wide  band  frequency 
when  it  comes  to  the  basic  functionality  of  an  antenna.  However,  narrow  band 
antennas  are  only  tuned  to  achieve  the  impedance  coupling  over  a  small  band  of 
frequency.  In  his  work,  “A  BRIEF  HISTORY  OF  UWB  ANTENNAS”,  Schantz 
(2004)  pointed  out  that,  ironically,  Oliver  Lodge  (1898)  has  mentioned  future  Ultra 
Wide  Band  antenna  structures  as  promising  designs  in  the  same  patent  which 
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introduced  the  concept  of  “syntony,”  or  the  idea  that  transmitter  and  a  receiver  should 
be  tuned  to  the  same  frequency  to  maximize  the  received  signal.  It  wasn’t  until  1939, 
when  Carter  introduced  his  conical  monopole  as  an  UWB  antenna  in  his  research  for  a 
wider  bandwidths  to  satisfy  the  demand  on  video  signals  increase.  Many  advances 
occurred  afterwards,  but  despite  the  excellent  performance  demonstrated,  other 
consideration  began  to  come  to  attention.  As  broadband  receivers  came  into  mass 
production,  the  demand  on  inexpensive  and  easy  to  manufacture  designs  increased. 
The  triangular  “bow-tie”  antenna  originally  proposed  by  Lodge  was  then  re-examined 
by  Brown  and  Woodward  in  1952  to  demonstrate  these  benefits. 

1.2.2.  Types  and  applications 

The  frequency  band  of  the  transmitted  signal  is  the  determining  factor  for  the 
GPR  type  of  application.  Higher  frequencies  are  often  needed  where  better  resolution 
is  required.  Nevertheless,  lower  frequencies  are  preferred  for  deeper  penetration  in 
soil  due  to  the  increasing  wave  attenuation  with  frequency  increase.  Ultra  Wide  Band 
(UWB)  GPR  systems  introduce  the  benefits  of  both  high  and  low  frequencies. 

In  general,  several  types  of  antennas  are  commonly  used  for  UWB  GPR 
systems,  such  as;  resistively  loaded  dipole,  bow-tie  antenna,  spiral  antenna  and  the 
most  widely  used  TEM  hom  antennas.  The  resistively  loaded  dipole  is  simple,  easy- 
to-design  and  linearly  polarized  antenna  structure.  However,  it  mostly  has  a  major 
disadvantage  of  poor  gain.  This  makes  other  planar  antennas,  like  Bow-tie  and  spiral 
structures,  a  more  attractive  choice. 

Bow-tie  and  spiral  antennas  are  mainly  used  in  impulse  GPR  applications 
because  of  their  non-dispersive  characteristics  and  their  relatively  high  gain. 
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However,  both  suffer  from  significant  impedance  mismatch  at  feed  point  when  used 


for  Air  Coupled  GPR  applications.  This  also  adds  the  ringing  effect  distorting  the 
time  domain  waveform.  Resistive  loading  is  usually  used  to  overcome  this 
drawback  but  this  technique  always  leads  to  deteriorating  gain  performance  over  the 
frequency  wide  band  (Lestari  et  al.,  2004).  Figure  3  and  Figure  4  show  the  typical 
structures  of  a  bow-tie  and  a  spiral  antenna  respectively. 


L 


L 


Figure  3:  An  example  of  Bow-Tie  antenna  structure 


Figure  4:  An  example  of  Spiral  antenna  structure 


1.2.3.  Horan  antennas 

TEM  Horn  antenna,  on  the  other  hand,  has  a  clear  advantage  over  planar 
antennas.  A  typical  TEM  Horn  antenna  structure  is  shown  in  Figure  5.  Its  narrow 
beam  width  naturally  guarantees  a  higher  directivity  gain  over  a  wider  frequency 
range.  The  main  design  challenge  is  to  achieve  a  good  impendence  match  in  order  to 
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minimize  feed  point  and  aperture  reflections.  The  work  in  Turk  et  al.  (2013)  is  an 
example  attempt  to  modify  the  TEM  horn  antenna  to  achieve  a  better  impedance 
matching.  The  TEM  horn-fed  ridged  horn  antenna  proposed,  demonstrates  a  high  gain 
up  to  lOdB  over  a  wide  frequency  band,  400MHz  to  16GHz,  with  input  reflections 
down  to  -lOdB.  However,  simulation  results  show  that  this  design  suffers  from 
relatively  high  reflections  over  frequencies  less  than  3GHz  and  greater  than  7GHz. 


A 


Figure  5:  An  example  of  a  Horn  antenna  structure 


1.3.  GPR  system  description  and  antenna  requirements 


We  have  been  developing  high  speed  impulse  ground  penetrating  radar  (Xianlei 
Xu  et  al.,  2012)  for  highway  pavement  and  bridge  deck  inspection.  This  GPR  system 
is  designed  to  achieve  high  spatial  resolution,  i.e.  1  cm  and  accurate  inspection  results 
while  scanning  at  regular  highway  speed,  i.e.  60  mph.  In  order  to  accomplish  high 
detection  resolution  while  reducing  power  radiation,  a  dual-band  dual-channel  UWB 
GPR  is  developed.  Figure  6  presents  the  GPR  system  diagram.  The  transmitter 
consists  of  an  UWB  pulse  generator  fed  by  a  tunable  1  kHz  to  10  MHz  oscillating 
signal  from  a  programmable  oscillator.  The  oscillating  signal  controls  the  UWB  pulse 
repetition  frequency  (PRF).  The  UWB  pulse  generator  provides  600ps-wide 
monocycle  pulses  to  both  base  band  and  3-4.5GHz  modulated  band.  In  the  base-band 
operation,  the  monocycle  pulse  is  directly  amplified  and  transmitted  while  at  the  high 
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frequency  band  the  monocycle  pulse  is  fed  to  a  mixer  for  frequency  modulation.  The 
selection  of  the  operating  band  is  remote-controlled  using  single-pole/dual-throw  RF 
switches.  The  receiver,  on  the  other  hand,  consists  of  the  receiving  antenna  flowed  by 
low  noise  amplifier  (LNA).  The  amplified  signal  is  either  directly  fed  to  the  analog- 
to-digital  (ADC),  or  frequency  down  converted  through  a  mixer,  depending  on  the 
operation  band.  The  digital  data  is  then  stored  in  a  memory  block  for  further 
processing.  To  improve  radar  survey  accuracy,  two  operating  channels  are 
implemented  for  parallel  inspections.  The  FPGA  coordinates  the  radar  operating  band 
selection,  channel  switching,  sampling  data  collection,  and  storage.  As  an  important 
functional  unit,  the  antenna  plays  a  significant  role  in  determining  GPR  performance. 
UWB  GPR  antenna  design  is  very  challenging  because  it  requires  exceptional 
impedance  matching  across  the  whole  ultra-wide  bandwidth. 


Programmable 

Oscillator 

UWB  Pulse 
Generator 

Povrtr 

Amplifier 


Figure  6:  Block  diagram  of  the  high-speed  UWB  dual-band  GPR  system  (Xianlei  Xu  et  al.,  2012) 
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CHAPTER  2:  ANTENNA  STRUCTURE 
2.1.  Antenna  sections 

Impedance  matching  is  one  of  the  most  critical  parts  in  UWB  antenna 
design.  Unfortunately,  it  is  difficult  to  achieve  an  exact  match  along  the  whole  wide 
band.  Basically,  there  are  two  main  points  needing  intensive  considerations,  one  is 
the  feeding  port  of  the  transmitter,  and  the  other  is  the  interface  between  antenna 
aperture  and  air.  For  the  air  coupled  GPR  system,  the  loading  impedance  of  the 
aperture  is  equal  to  the  characteristic  impedance  of  air  which  can  be  calculated 
using  equation  (4).  However,  given  that  the  air  electric  conductivity  a  is  very  small 
in  value  and  both  pr  and  sr  can  be  approximated  to  1,  the  characteristic  impedance 
can  then  be  approximated  to: 

Z0  =  =  VmoAo  ~  377  n 

On  the  other  hand,  the  impedance  of  the  feeding  coaxial  cable  is  50  Ohms.  A 
good  antenna  design  should  resolve  this  mismatch  by  gradually  converting  the 
impedance  to  minimize  reflections  and  therefore  minimize  the  Voltage  Standing 
Wave  Ratio  (VSWR).  In  return  this  leads  to  maximizing  radiation  efficiency.  Figure 
7  and  Figure  8  show  that  our  new  design  consists  basically  of  three  sections;  feed 
line,  Vivaldi  shaped  section  and  rounded  shaped  aperture.  Essentially  our  proposed 
design  combines  the  benefits  of  both  Vivaldi  and  horn  antennas. 
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Figure  7:  Top  view  of  Vivaldi  Shaped  Horn  antenna 


Figure  8:  Side  view  of  Vivaldi  Shaped  Horn  antenna 


2.2.  Analytical  model 

Both  feed  line  and  the  linear  Vivaldi  section  can  be  modeled  as  combination 
of  N  parallel  plate  transmission  line  segments.  Each  segment  consists  of  two  metal 
plates,  with  varying  width  and  space,  separated  by  a  dielectric. 
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Figure  9:  Feed  line  as  parallel-plate  transmission  line 

Starting  with  the  feed  line,  the  relationship  between  input  and  output 
impedances  is  controlled  by  equation  (5). 


y  _  y  ZOUt  +  yZ0tan(/?Z) 
in~  °Z0+jZouttan(pi) 

Where: 


(5) 


Zin  input  impedance  of  the  parallel  plate  wave  guide 
Zout  output  impedances  of  the  parallel  plate  wave  guide 
/  The  length  of  the  feed  line. 

Zq  The  characteristic  impedance 

[j  The  wave  number  of  the  parallel  plate 


The  characteristic  impedance  Zo  and  the  wave  number  [>  of  the  transmission 
line  are  defined  by  equations  (6)  and  (7)  where  d  is  the  dielectric  height  that 
separates  the  two  parallel  plates  and  w  is  the  width  of  the  feed  line. 


(6) 


/?  = 


(7) 
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In  order  to  get  a  better  understanding  of  the  parallel  plate  structure,  the 
MATLAB  script,  presented  in  APPENDIX  section,  is  developed  to  get  the  3D  plot 
of  the  variations  of  input  impedance  versus  the  change  in  both  frequency  and  the 
length  of  our  feed  point.  Given  an  output  impedance  of  50  ohms  and  a  dielectric  of 
air,  Figure  10  and  Figure  11  show  the  real  and  imaginary  parts  of  the  input 
impedance  respectively. 


Figure  10:  Real  part  of  input  impedance  versus  length  and  frequency 


20  0 


Figure  11:  Imaginary  part  of  input  impedance  versus  length  and  frequency 
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Tabulated  results  are  shown  in  Table  1,  where  the  input  impedance  is 


reported  for  each  length  value  along  its  corresponding  frequency  points.  The 
analysis  shows  that  when  the  feed  line  has  dimensional  size  of  LO  =  6mm,  and  dO  = 
3mm  and  WO  =  12mm,  the  feed  line  demonstrates  minimal  impedance  variation 
across  the  wide  frequency  band. 


Table  1  -  Input  impedance  for  different  length  values  over  frequency  band  0.5GHz  to  6.0GHz 


L0=  6  mm 

L0=  9  mm 

L0=  12  mm 

L0=  18  mm 

f=  0.5GHz 

50.14 +  4.250i 

50.32  +  6.380i 

50.57  +  8.510i 

51.30+  12.78i 

f=  1.0GHz 

50.57 +  8.510i 

51.30+  12.78i 

52.33  +  17.06i 

55.40  +  25.66i 

f=  1.5GHz 

51.30+  12.78i 

52.96+  19.21i 

55.40  +  25.66i 

63.01  +  38.58i 

f=  2.0GHz 

52.33+  17.06i 

55.40  +  25.66i 

60.02  +  34.29i 

75.40  +  50.90i 

f=  2.5GHz 

53.69  +  21.36i 

58.70 +  32.14i 

66.53 +  42.8H 

94.43  +  60.72i 

./=  3.0GHz 

55.40 +  25.66i 

63.01  +  38.58i 

75.40  +  50.90i 

121.48  +  63.23i 

f=  3.5GHz 

57.50  +  29.98i 

68.50  +  44.89i 

87.23  +  57.94i 

153.4  +  49.82i 

f=  4.0GHz 

60.02  +  34.29i 

75.40  +  50.90i 

102.55  +  62.72i 

175.7  +  14.65i 

f=  4.5GHz 

63.01  +  38.58i 

83.96  +  56.34i 

121.48  +  63.23i 

170.46  -  28.93i 

/=  5.0GHz 

66.53 +  42.8H 

94.43  +  60.72i 

142.84  +  56.65i 

142.46  -  56.84i 

f=  5.5GHz 

70.63  +  46.94i 

106.96  +  63.35i 

162.97  +  40.38i 

lll.26-63.66i 

./=  6.0GHz 

75.40 +  50.90i 

121.48  +  63.23i 

175.7  +  14.65i 

87.00  -  57.83i 

The  other  perspective  that  should  be  of  interest  is  the  propagation  modes  that 
travel  through  the  parallel  plate  structure.  Given  that  the  desired  antenna  will 
operate  at  a  wide  range  of  frequencies,  it  is  important  to  understand  the  modes 
traveling  through  the  antenna  at  various  frequencies.  The  TEM  propagation  mode  is 
the  dominant  as  it  has  a  cutoff  frequency  of  zero.  However,  TE  and  TM  modes  also 
exist  but  with  higher  cutoff  frequencies.  Equation  (8)  shows  that  for  each 
propagation  mode  n,  the  cutoff  frequency  fcn  is  inversely  proportional  to  the 
dielectric  thickness  d ,  the  square  root  of  magnetic  permeability  p  and  the  square 
root  of  the  electrical  permittivity  s. 
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fm  2  diyfjli)  Vn  0’1,2’- 


(8) 


At  the  feed  point,  d  is  at  its  minimum  value,  and  equals  to  dO  where  dO  = 
3mm.  Equation  (8)  can  directly  be  used  to  calculate/,,.  However,  the  local  dielectric 
thickness  dA  at  the  aperture,  the  end  of  the  Vivaldi  section,  is  calculated  in  equation 
(9)  where  Vieng]lt  the  length  of  the  Vivaldi  section  assumed  to  be  180mm  and 
D_angle  is  set  to  5.5°. 

dA  =  d0  +  Vlength  x  sin(D_angle )  (9) 

By  substituting  in  equation  (8),  the  cutoff  frequency  at  aperture  fcnA  is 
calculated  by  equation  (10). 


fcnA 


n 


(d0 


+  V, 


length 


_ (J_) 

x  sin  (Dan„le))  ^ 


(10) 


Table  2  shows  the  calculated  values  of  each  cutoff  frequency  at  both  ends  of 
the  antenna. 


Table  2:  Cutoff  frequencies  of  each  propagation  mode  at  feed  point  and  aperture 


Fc 

Feed 

aperture 

Fcl 

49.97  GHz 

3.99  GHz 

Fc2 

99.93  GHz 

7.97  GHz 

Fc3 

149.90  GHz 

11.96  GHz 

Fc4 

199.86  GHz 

15.93  GHz 

Fc5 

249.83  GHz 

19.92  GHz 
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Figure  12:  Cutoff  frequencies  of  each  propagation  mode 


For  the  Vivaldi  section,  it  can  be  represented  by  staircase  model  shown  in 
Figure  13,  which  consists  of  N  segment.  Each  segment  can  be  assumed  homogenous 
when  the  segment  length  is  small  in  comparison  with  signal  wavelength. 


Figure  13:  Staircase  model  of  the  Vivaldi  section 

Similarly,  the  input  impedance  Zin  of  each  stage  can  then  be  calculated 
using  equation  (11),  where  pi,  ZOi  and  li  are  the  wave  number,  the  characteristic 
impedance  and  the  length  of  the  i’th  segment  respectively.  Zini+1  is  the  input 
impedance  of  the  i+1  stage  that  loads  the  i’th  segment.  The  point  where  i  =  0, 
represent  the  touching  point  of  the  feed  line  and  Vivaldi  sections.  Where  i  =  N 
represent  the  other  end  of  the  Vivaldi  section  where  it  touches  the  arc  section. 
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An  analytical  model  developed  in  MATLAB  is  created  based  on  Equations 
(11)  and  (12)  to  find  suitable  values  of  design  angels;  D_angle  and  W_angle  shown 
in  Figure  7  and  Figure  8  respectively.  The  code  is  available  in  the  APPENDIX 
section.  In  the  model,  the  loading  impedance  of  the  last  stage  Zl  is  equal  to  free 
space  characteristic  impedance  (377  Ohms).  The  Nth  segment  of  the  Vivaldi  section 
is  chosen  to  interface  directly  with  air,  dropping  out  the  arc  section  for  simplicity. 
Equation  (12)  is  then  used  to  plot  SI  1  behavior  at  the  feed  point. 

D_angle  is  swept  from  2.75  to  5.5°  with  a  step  of  0.25°  and  W_angle  is 
swept  from  13°  to  20°  with  a  step  of  0.5°.  Figure  14  shows  the  Sll  results  at  the 
start  and  end  points  of  each  sweep  range  for  simplicity.  It  shows  that  when  D_angle 
at  its  largest  value  5.5°  and  W_angle  at  its  smallest  value  13°,  promising  Sll  results 
across  the  wide  frequency  range  is  achieved. 
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Figure  14:  Sll  results  for  a  preliminary  MATLAB  model  of  feed  point  and  Vivaldi  sections 
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CHAPTER  3:  FIELD  SOLVER  SIMULATION 


3.1.  Model  construction 


Starting  by  building  the  structure  in  SolidWorks,  the  basic  design  parameters 
are  initially  plugged  to  create  the  mechanical  model.  This  helped  later  in 
manufacturing  this  antenna  structure  to  perform  lab  measurements.  The  second  step 
is  to  import  the  model  into  HFSS  for  design  parameter  optimization.  Feed  port  is 
added  and  the  simulation  setup  is  prepared  to  obtain  the  far  field  simulation  results. 
HFSS  simulation  results  guides  us  through  tuning  up  different  parts  of  the  antenna 
to  obtain  the  optimum  performance  of  this  structure. 

Figure  15  illustrates  that  the  exact  model  of  the  coaxial  signal  probe  with  0.5 
inch  flange  is  first  imported  in  SolidWorks  to  remove  the  outer  threading  and  to 
shorten  the  probe  length.  The  threading  removal  helps  later  when  the  model  is 
imported  in  HFSS  for  simulation.  The  thread  is  a  complex  feature  that  requires  more 
convergence  timing  in  electromagnetic  field  simulation. 


Figure  15:  Coaxial  feed  modeling 


After  importing  the  antenna  structure  in  HFSS,  we  place  a  vacuum  box  around. 
The  walls  of  the  box  are  selected  as  radiation  boundaries  for  intercepting  the  far  field 
radiation.  The  size  of  the  box  has  to  be  greater  than  the  longest  feature  in  the  antenna 
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structure,  in  our  case  it’s  the  antenna  total  length  of  300mm.  This  is  to  satisfy  the  far 


field  radiation  conditions. 


Figure  16:  HFSS  3D  antenna  model 


The  excitation  is  then  assigned  to  face  of  the  dielectric  material  of  the  coaxial 
probe.  The  input  impedance  is  set  to  50mm.  This  matches  the  standard  connector 
impedance  that  connects  our  antenna  to  the  impulse  generator  when  building  up  the 
GPR  system. 


Figure  17:  Assigning  excitations  in  HFSS 


3.2.  Parameters  optimizations 

One  of  the  critical  parameters  of  this  antenna  design  is  the  length  of  the 
Vivaldi  section.  Parametric  simulations  are  performed  to  evaluate  the  effects  on 
Vivaldi  section  length  on  Sll  across  the  wide  band  frequency.  To  perform  a 
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frequency  sweep  in  HFSS,  we  need  to  “Add  Solution  Setup”  from  “Analysis  Setup” 
menu  and  then  we  select  “Add  Frequency  sweep”.  Figure  18  shows  the  sweep  setup 
of  a  frequency  range  of  100MHz  to  6GHz  with  a  linear  step  lMFIz. 


Figure  18:  Frequency  Sweep  Setup  in  HFSS 


S 1 1  results  can  be  displayed  through  creating  a  report  using  “Create  Model 
Solution  Data  Report”  from  “Results”  menu  and  then  choosing  “Rectangular  Plot”. 
The  dialog  in  Figure  19  is  used  to  plot  Sll  in  dB  units  over  the  frequency  sweep 
that  is  called  “Sweep  1”. 


Figure  19:  Create  a  rectangular  plot  for  Sll  parameters  in  HFSS 
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Figure  20  illustrates  the  simulation  results  when  the  length  is  linearly 
changed  starting  from  90mm  to  180mm  with  a  step  size  of  30mm.  Sll  is  plotted  for 
each  length  value  correspondingly.  The  results  show  that  180mm  length  results  in 
steady  Sll  curve  below  -lOdB  across  the  whole  frequency  band. 


Figure  20:  Sll  vs.  Frequency  for  different  antenna  lengths 

The  shape  and  length  of  arc  section  also  stands  as  a  critical  design  variable. 
Given  that  the  lowest  frequency  is  600MFIz,  the  full  wave  length  is  500mm.  The 
total  length  of  the  antenna  is  chosen  to  be  about  half  the  wave  length.  Hence,  the  arc 
length  and  the  length  of  the  Vivaldi  section,  180mm,  should  add  up  to  250mm  to 
match  the  half-wavelength  of  the  lowest  end  frequency  of  the  frequency  band. 
Simulation  result  confirms  that  the  initial  arc  length  value  of  70mm  is  proven  to  be 
optimum. 
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Sll  (dB) 


Enhanced  performance  is  achieved  by  rounding  the  edges  at  the  feed  point. 
Figure  21a  and  Figure  21b  show  the  traditional  sharp  edges  versus  the  edge 
rounding  at  feed  point  respectively.  Results  in  Figure  22  demonstrates  that  the 
structure  in  Figure  21b  has  lower  Sll  results,  dashed  line,  over  the  solid  line  that 
represents  simulation  results  of  the  traditional  sharp  cornered  design. 


Figure  21:  Antenna  design,  feed  line  tune-ups 


Figure  22:  Sll  simulation  results  for  antenna  structures  in  Figure  21a  (solid  line)  vs.  Figure  21b 

(dashed  line) 
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Rounding  the  flare  edges  of  the  antenna  also  demonstrates  an  enhancement 
to  the  final  Sll  results.  Figure  25  shows  that  rounding  the  sharp  corners,  as  in 
Figure  24,  has  lower  Sll  results  (dashed  line)  versus  the  sharp  edges,  as  in  Figure 
23,  over  the  low  frequency  range  600MHz  to  1GHz.  Although  the  upper  end,  5GHz 
to  6GHz,  has  higher  Sll  results,  it  is  still  well  below  -lOdB  which  is  acceptable.  In 
conclusion  the  structure  in  Figure  24  achieved  better  performance. 


Figure  23:  Antenna  design  tune-ups.  Feed  line  modifications 


Figure  24:  Sll  simulation  results  for  antenna  structure  in  Figure  23a 
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(8P)  US 


Figure  25:  Sll  simulation  results  for  antenna  structure  in  Figure  23b 


3.3.  Simulation  results 


Figure  25  shows  the  radiation  pattern  for  the  final  antenna  design  at  the 


midrange  2. 8GHz  frequency  and  cp  =  0  degrees.  Figure  12b  shows  the  3D  polar  plot 


at  the  same  frequency  with  a  high  peak  gain  of  10.7dB  and  a  radiation  efficiency 


reaching  up  to  98%. 
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Figure  26:  Radiation  pattern  in  the  x-y  plan  (2.8  GHz) 
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Figure  27:  3D  polar  plot  (2.8  GHz) 
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CHAPTER  4:  EXPERIMENTAL  RESULTS 


4.1.  Prototype  manufacturing 

After  updating  the  SolidWorks  mechanical  model  with  the  optimum  design 
parameters  verified  by  HFSS,  a  first  prototype  is  then  manufactured.  Figure  28 
shows  the  exact  dimensions  of  one  of  the  antenna  flattened  copper  plates.  Copper 
sheet  with  0.5mm  thickness  is  cut  to  the  pattern  below.  The  two  plates  are  then 
bended  at  the  end  of  the  Vivaldi  section  to  form  an  arc  with  a  120mm  length  as  a 
part  of  a  complete  circle  with  45mm  radius.  Rounding  the  edges  of  both  the  feed 
line  and  aperture  was  also  performed  based  of  the  SolidWorks  mode.  The  coaxial 
connector  is  welded  between  the  two  plates,  maintaining  the  3mm  gap  in  the 
antenna  design  and  providing  some  robustness  at  the  feed  line. 


12.21 


„ _ 10.29 

(.47 


Figure  28:  Antenna’s  right  copper  plate 


To  enhance  the  robustness  of  the  structure,  the  antenna  is  surrounded  by  a 
material  with  dielectric  constant  equals  to  approximately  1,  in  order  not  to  affect 
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performance.  Figure  29  shows  the  final  shape  of  the  prototype  manufactured  in 
UVM  laboratory.  The  structure  is  then  placed  in  a  regular  cardboard  box  as  shown 
in  Figure  31  for  field  measurements. 


Figure  29:  First  prototype  antenna 

4.2.  Lab  results 

Sll  measurement  is  then  performed  using  HP  8753D  network  analyzer.  The 
magnitude  of  the  reflection  coefficient  is  measured  over  a  frequency  span  from 
100MHz  to  6GHz  with  a  20MHz  step.  Figure  30  shows  the  plot  is  perfectly 
matching  the  simulation  results  swept  over  the  same  frequency  range  with  a  1MHz 
step. 

The  test  setup  in  Figure  31  uses  a  pair  of  transmitting  and  receiving  antennas 
to  measure  the  time  domain  response.  The  transmitted  signal  is  a  Gaussian 
monocycle  pulse;  that  has  pulse  width  of  approximately  600ps  and  amplitude  of 
approximately  1.8V.  Our  second  prototype,  the  receiver,  is  placed  exactly  1  meter 
away  from  the  transmitter  and  connected  to  the  spectrum  analyzer  to  capture  the 
waveform  of  the  signal  received. 
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Figure  30:  Sll  measurements  vs.  simulation  results  of  the  prototype  antenna 


Figure  31:  Time  domain  waveform  test  setup 
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The  same  experiment  is  repeated  using  a  pair  of  A.H.  Systems’  double  ridge 
guide  hom  antennas  model  S AS -571  for  performance  comparison.  Figure  31  shows 
that  using  our  antennas,  the  received  pulse  features  smaller  ringing  effect  than  using 
the  double  ridge  guide  hom  whose  pulse  signal  received  is  shown  in  Figure  30.  Many 
factors  can  contribute  to  this  better  performance;  including  linear  phase  and  steady 
high  gain  across  the  wide  frequency  band.  However,  we  believe  that  the  major  reason 
is  the  good  matching  achieved  at  lower  frequencies,  less  than  1  GHz. 


Figure  32:  Transmitted  impulse 
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Figure  33:  Received  signal  from  A.H.  Systems’  Double  Ridge  Guide  Horn  Antennas  (model  SAS- 

571) 


Figure  34:  Received  signal  from  our  antenna 
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4.3.  Field  results 


In  order  to  measure  the  far  field  radiation  pattern,  the  test  setup  in  Figure  35 
uses  an  isotropic  antenna  as  a  receiver  placed  lm  apart  from  our  transmitter.  The 
transmitting  antenna  is  placed  over  a  rotating  table  that  allows  us  to  measure 
radiating  power  over  E-plan  angles  from  0°  to  180°  with  a  step  of  10°.  The  open 
area  test  site  is  equipped  with  buried  grounded  mesh  that  acts  as  a  perfect  ground 
plane.  However,  for  more  accurate  readings,  the  presented  results  reflect  the  same 
experiment  conducted  in  an  electromagnetic  compatibility  (EMC)  lab.  The 
laboratory  is  a  concealed  chamber  that  has  electrostatic  discharge,  magnetic  and 
radiation  immunity  capabilities.  The  input  signals  are  fed  using  standard  50  Ohms 
cables  through  chamber  outlets. 

For  a  single  frequency,  the  reading  of  the  received  RF  voltage  VRF  is 
captured  in  (Volts)  and  converted  to  (dB  V)  using  equation  (7). 

V_RF  (dB  V)  =  20  log io  [V_RF  (V)]  (7) 

The  losses  due  to  cables  and  connectors  are  determined  by  connecting  the 
cable  under  test  between  the  signal  generator  and  the  spectrum  analyzer  at  the 
frequency  of  interest,  and  examine  how  much  attenuation  the  signal  suffered. 
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Figure  35:  Radiation  pattern  test  setup 


On  the  other  hand,  the  antenna  calibration  process  starts  with  replacing  our 
antenna  with  A.H.  Systems’  (SAS-571)  Double  Ridge  Guide  Horn  Antenna  and 
measure  the  power  received  at  the  same  frequency  points.  The  readings  are  then 
compared  to  the  tabulated  expected  readings  documented  in  the  SAS-571  datasheet. 
The  difference  is  then  subtracted  from  our  antenna  measurements  to  acquire  the 
calibrated  radiation  pattern.  Figure  36  and  Figure  37  show  that  our  design  has 
demonstrated  a  very  good  match  relative  to  simulation  results  over  the  frequencies 
of  1.9  GHz  and  2.8  GHz.  The  solid  line  on  both  figures  represents  the  simulated 
radiation  pattern  over  E-plan  angles  from  0°  to  360°  with  a  step  of  1°. 
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Figure  36:  Measured  (dashed)  versus  simulation  (solid)  radiation  pattern  at  1.9GHz 


Figure  37:  Measured  (dashed)  versus  simulation  (solid)  radiation  pattern  at  2.8GHz 
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CHAPTER  5:  CONCLUSION 


5.1.  Goals  achieved 

The  design  procedure  from  initial  concept  to  manufacturing  the  first 
prototype  has  been  demonstrated.  Both  targets  of  compact  design  dimensions  and 
wideband  operation  have  been  achieved.  Throughout  the  design  process  we  have 
identified  the  critical  physical  features  to  tune  in  order  to  achieve  the  optimum 
results.  The  antenna  characteristics  have  been  simulated  and  measured.  Its  good 
performance  has  been  verified  in  several  test  setups.  Our  antenna  is  suitable  for  the 
detection  of  subsurface  structure  with  good  resolution. 

During  the  course  of  this  thesis,  a  design  methodology  has  been  developed 
to  test  any  horn  antenna  structure  by  first  changing  the  contributing  dimensions  in 
the  analytical  MATLAB  model  to  get  a  proof  of  concept  if  these  dimensions  will 
have  enough  potential  to  proceed  in  the  next  design  step.  The  second  step  is  to 
match  the  SolidWorks  model  to  the  desired  dimensions  and  finally  import  the 
structure  in  HFSS  for  field  simulation.  Design  iterations  will  guarantee  the  design 
readiness  for  production. 

5.2.  Future  work 

The  modeling  technique  used,  opens  the  door  to  many  possibilities.  The  first 
and  foremost  near  objective  is  to  integrate  the  antenna  developed  with  the  GPR 
system  for  production.  All  the  field  testing  is  already  performed  using  the  blocks 
from  the  system  proposed,  namely  the  UWB  impulse  generator,  the  LNA  and  the 
mixer  for  high  frequency  band  operation.  However,  good  packaging  and  handling 
the  system  reliability  is  still  requires  a  follow  up. 
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This  thesis  provides  a  bundle  of  Swiss-knife  set  of  software  tools  that  can  be 
helpful  to  perform  task  specific  operation.  Many  enhancements  can  be  added  to 
each  software  module  for  shorter  simulation  time,  better  accuracy  or  a  more  user- 
friendly  interface.  The  set  of  MATLAB  scripts  included  represents  preliminary 
analytical  models  that  can  be  enhanced  to  improve  result  accuracy.  More  models 
can  also  be  added  to  this  set  to  include  2D  antennas  in  addition  to  the  current  only 
supported  Horn  antennas.  A  graphical  user  interface  can  also  be  developed  to  point 
out  the  different  input  features  to  the  model,  so  it  becomes  more  user-friendly. 

One  of  the  most  challenges  in  this  work  is  the  communication  link  between 
SolidWorks  and  HFSS.  The  increase  usage  of  3D  printers  and  the  powerful 
capabilities  of  SolidWorks  as  solid-modeling  CAD  software,  demands  an  easy-to- 
use  link  to  powerful  software  in  electromagnetic  simulation  similar  to  HFSS.  A 
program  that  runs  both  tools  in  the  background  and  provide  the  result  in  GUI  format 
is  definitely  a  nice  to  have. 
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APPENDIX 


The  following  MATLAB  module  is  used  to  calculate  input  impedance  for  a 
parallel-plate  transmission  line  with  a  length  l  and  output  impedance  of  50  Ohms. 

%  Start  a  clean  session 

clear; 

close  all; 

clc; 

zout  =  50;  %  output  impedance 

d  =  0.003;  %  fixed  dielectric  hight 

w  =  0.012;  %  fixed  dielectric  width 

meu  =  4*pi*10A-7;  %  meuO  of  free  space 

eps  =  8.854187817  *  10A-12;  %  epsO  of  free  space 

zO  =  (d/w)  *  sqrt(meu/eps);  %  characteristic  impedance  of  free  space 

%  Loop  on  frequency  range  starting  at  1MHz  and  ends  at  6GHz, 

%  with  a  step  of  10MHz 

j=0; 

for  f  =  1*10A6  :  10A7  :  6*10A9 

j=j+i; 

freq(j)  =  f/10A6; 
k  =  0; 

%  Loop  on  transmission  line  length,  starting  from  6mm  to  18mm, 
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%  with  a  step  of  1mm 
for  1  =  0.006  :  0.001  :  0.018 
k=k+l; 

beta  =  2*pi*f*sqrt(meu*eps);  %  Calculate  beta 

%  Calculate  input  impedance 

Zin  =  zO  *((zout+li*zO*tan(beta*l))  /  ((zO+li*zout*tan(beta*l)))); 
length(k)  =  1*10A3; 

RZin(j,k)  =  real(Zin); 

IZin(j,k)  =  imag(Zin); 
end 
end 

%  Plot  real  part  of  input  impedance 
me  she  (length, freq,RZin); 
xlabel  ('Length  (mm)'); 
ylabel  ('Frequency  (MHz)'); 
zlabel  ('Real  Impedence  (Ohms)'); 

%  Open  a  new  figure 
figure; 

%  Plot  imaginary  part  of  input  impedance 
meshc  (length, freq,IZin); 
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xlabel  ('Length  (mm)'); 

ylabel  ('Frequency  (MHz)'); 

zlabel  ('Imaginary  Impedence  (Ohms)'); 
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The  following  is  the  MATLAB  code  is  used  to  calculate  the  final  input 
impedance  when  combining  the  feed  line  with  the  Vivaldi  section  designed. 

%  Start  a  clean  session 

clear; 

close  all; 

clc; 


dielectric  hight 
dielectric  width 

%  feedline  dielectric  meu  =  meuO 
%  feedline  dielectric  eps  =  epsO 
%  speed  of  light  in  free  space 
%  ~  377  Ohms 


dO  =  0.003;  %  feedline 

wO  =  0.012;  %  feedline 

meu  =  4*pi*10A-7; 
eps  =  8.8541 878 17*  10^- 12; 
c0  =  3*10A8; 

Z0_air  =  sqrt(meu/eps); 

D_angle  =  5.52*pi/180; 
W_angle  =  12.95*pi/180; 


%  Vivaldi  section  D_angle 
%  Vivaldi  section  W_angle 


AnLength  =  0.1 80;  %  Assuming  a  lenght  of  1 80mm  for  Vivaldi  section 

j  =  0; 

%  Loop  on  frequency  range  starting  at  1MHz  and  ends  at  6GHz, 

%  with  a  step  of  10MHz 
for  f  =  1*10A7  :  1*10A7  :  6*10A9 
lamda  =  cO/f;  %  Calculate  lamda 
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1  =  lamda/1000; 


beta  =  2*pi*f*sqrt(meu*eps);  %  Calculate  beta 

%  Construct  a  vector  with  antenna  Vivaldi  section  length 
V_length  =  AnLength:-l:0; 

%  Construct  a  vector  with  antenna  Vivaldi  section  local  dielectric  thickness 
V_d  =  dO  +  V_length.*sin(D_angle).*2; 

%  Construct  a  vector  with  antenna  Vivaldi  section  local  dielectric  width 
V_w  =  wO  +  V_length.*sin(W_angle).*2; 

zout  =  Z0_air;  %  Initialize  loading  impedance  with  air  intrinsic  impedance 

%  Loop  over  each  stage  of  the  staircase  model 
for  k  =  1  :  length(V_length) 
if  (k==length(V_length)) 

%  If  it’s  the  last  stage  towards  the  feed  line: 
d  =  (dO+V_d(k))/2;  %  get  the  average  of  the  d  value 

w  =  (wO+V_w(k))/2;  %  get  the  average  of  the  w  value 

1  =  V_length(k); 
else 

%  If  it’s  any  other  stage  starting  from  the  aperture  towards  the  feed  line: 
d  =  (V_d(k+l)+V_d(k))/2;  %  get  the  average  of  the  d  value 
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w  =  (V_w(k+l)+V_w(k))/2; 


%  get  the  average  of  the  w  value 


end 

%  Calculate  intrinsic  impedance  for  each  stage 
zO  =  (d/w)  *  sqrt(meu/eps); 

%  Calculate  input  impedance  for  each  stage 

Zin(k)  =  zO  *((zout+li*zO*tan(beta*l))/((zO+li*zout*tan(beta*l)))); 

%  Assign  current  input  impedance  as  the  output  impedance  for  the  next  iteration 
zout  =  Zin(k); 
end 

%  This  is  the  input  impedance  of  the  last  stage  in  Vivaldi  section 
Zin_vivaldi  =  zout; 

%  Local  parameters  for  the  feed  line 
feed_d  =  0.003; 
feed_w  =  0.012; 
feed_l  =  0.012; 

feed_z0  =  (feed_d/feed_w)  *  sqrt(meu/eps); 

%  Calculate  input  impedance  for  the  feed  line 
Zin  =  feed_z0  *  ((Zin_vivaldi+li*feed_zO*tan(beta*feed_l))  / 
((feed_zO+li*Zin_vivaldi*tan(beta*feed_l)))); 
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j=j  +  i; 


freq(j)  =  f; 

%  Final  input  impedance  for  this  frequency  node 
ABS_Zin(j)  =  abs(Zin); 

ANG_Zin(j)  =  angle(Zin)*  180/pi; 

%  Final  S 1 1  value  for  this  frequency  node 
SI  10)  =  20*logl0(abs((Zin-50)/(Zin+50))); 
end 

%  plot  input  impedance  in  absolute  and  angle 

[AX,  HI,  H2]  =  plotyy  (freq,ABS_Zin,freq,ANG_Zin); 

grid  on; 

%  set  axes  limits  and  label 
axes  (AX(1)) 
ylim([-450  550]) 
ylabel  ('Zin  Aboslute'); 

axes  (AX(2)) 
ylim([-90  90]) 
ylabel  ('Zin  Angle'); 
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xlabel  ('Frequency'); 

title(['Vivaldi  length  is  n u m 2str  ( AnLength*  1 000) ,  'mm']) 

%  Open  a  new  Figure 
figure; 

%  Plot  S 1 1  results 
plot  (freq,Sll,'r'); 
ylabel  ('Sll  (dB)'); 
xlabel  ('Frequency  (Hz)'); 
grid  on; 
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The  following  MATLAB  code  is  used  to  calculate  the  cutoff  frequency  for 
propagation  modes  in  a  parallel-plate  transmission  line. 

%  Start  a  clean  session 
clear;  close  all;  clc; 

dO  =  0.003;  %  fixed  dielectric  hight 

w  =  0.012;  %  fixed  dielectric  width 

meu  =  4*pi*10A-7;  %  meuO  of  free  space 

eps  =  8.854187817  *  10A-12;  %  epsO  of  free  space 

V_length  =  0.180;  %  Assuming  a  lenght  of  180mm  for  Vivaldi  section 

D_angle_variable  =  5.52;  %  Vivaldi  section  D_angle 

D_angle  =  D_angle_variable*pi/180;  %  D_angle  in  radians 

%  aprature  dielectric  thickness 
dl  =  dO  +  V_length*sin(D_angle)*2; 

%  loop  over  propagation  modes  "n"  from  1  to  5 
for  n  =  1  :  5 

%  cutoff  frequency  at  feed  line 
Fon_0(n)  =  n/(2*d0*sqrt(meu*eps)*10A9); 
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%  cutoff  frequency  at  aprature 


Fon_l(n)  =  n/(2*dl*sqrt(meu*eps)*10A9); 


end 


%  round  values  to  the  second  dicimal  point 
finall  =  (round  (Fon_0*  100)/ 100) 
fina!2  =  (round(Fon_l*100)/100) 
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